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1.  Introduction 


1.1  Background 

In  the  late  1980's  it  was  determined  that  small  relativistic  electron  beam  linear  accelerators  (linac's)  had 
been  reduced  in  size  to  the  point  where  they  could  be  feasibly  flown  onboard  spacecraft  and  sub-orbital 
rockets  [Banks  etal,  1987;  Banks  et  al.  1990].  Ionospheric  modification  by  beams  of  this  energy  can 
create  significant  ionization  and  steep  potential  gradients  resulting  from  the  interaction  of  the  relativistic 
(MeV)  beams  with  the  surrounding  environment.  These  beams  should  be  capable  of  achieving  deep 
penetration  into  the  lower  ionosphere  and  even  into  the  upper  stratosphere  when  fired  directly 
downwards  along  the  Earth's  magnetic  field  from  an  orbital  or  sub-orbital  spacecraft. 

Earlier,  the  Air  Force  sponsored  preliminary  research  with  the  University  of  Michigan  to  assess  the 
practical  aspects  of  relativistic  beam  propagation  including  the  extent  a  relativistic  electron  beam  was 
scattered  laterally  by  collisions  with  the  neutral  atmosphere  as  it  propagates  downwards  into  the 
atmosphere  from  a  space-based  platform.  It  was  found  that  the  beam  is  highly  collimated  by  the 
magnetic  field,  and  therefore  well  suited  for  ionospheric  modification-type  experiments.  This  study 
discussed  a  number  of  scientific  investigations  that  could  be  done  with  relativistic  electron  beams.  They 
include:  (1)  Study  of  the  energetic  electron  population  in  the  magnetosphere,  (2)  Modification  of  the 
atmospheric  electric  potential  structure,  and  (3)  Modification  of  chemical  reaction  paths  in  the  middle 
atmosphere  stimulated  by  naturally  precipitating  relativistic  electrons.  In  a  later  study  Neubert  et  al. 
[1996],  also  showed  that  beam  confinement  by  the  external  axial  magnetic  field  significantly  reduces  the 
amount  of  radial  expansion  due  to  small-angle  scattering  and  inelastic  collision  losses,  constraining  the 
beam’s  radius  to  be  less  than  its  gyroradius  at  a  90®  pitch  angle.  These  initial  results  provided  sufficient 
encouragement  to  scientists  at  the  Air  Force  Research  Laboratory  (AFRL)  to  consider  that  long-range 
propagation  of  MeV  electrons  beams  in  the  Earth’s  upper  atmosphere  is  theoretically  feasible. 
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The  Spacebome  Linear  Accelerator  System  (SLESfAC)  is  a  program  designed  to  demonstrate  the  viability 
of  injecting  relativistic  (E  -  MeV)  electron  beams  into  the  earth's  upper  atmosphere/ionosphere. 
Motivating  forces  behind  this  mission  include  [Jost,  1993]: 

•  Upper  atmosphere/ionosphere/magnetosphere  research 

•  Orbital  debris  detection  and  mitigation 

•  Atmospheric  ozone  depletion  monitoring  and  mitigation 

•  RF  scattering  screens 

•  Department  of  Energy  directed  beam  considerations 

The  beam  generator  developed  for  the  Air  Force  in  this  program  has  the  following  characteristics  that  are 
used  here  for  practical  considerations: 

•  beam  energy:  5  MeV 

•  beam  current:  100  mA 

•  pulse  length:  4  ^is 

•  duty  cycle:  0.001 

1.2  Goals 

The  overall  goal  of  the  work  of  this  grant  was  to  develop  computational  tools  and  use  them  to  better 
assess  relativistic  electron  beam  launch  and  propagation  from  space  platforms  and  its  interaction  with  the 
space  environment  and  atmosphere.  To  accomplish  this  goal,  several  computational  tools  were  developed 
and  then  applied  to  this  problem.  These  include: 

1.  The  paraxial  envelope  model  of  beam  propagation:  A  set  of  one-dimensional  equations  to 
calculate  the  beam  radius  and  energy  was  applied  to  propagation  in  the  atmosphere,  ionosphere, 
and  magnetosphere.  The  model  was  also  used  to  support  an  analysis  of  various  instability 
growth  rates. 

2.  Relativistic  Beam-Atmosphere  Interaction  (RelBAI)  model:  A  model  based  on  the  two- 
stream  approach,  with  extension  to  relativistic  energies  using  the  paraxial  envelope  equations 


was  developed  to  model  energy  degradation  and  optical  emissions  resulting  from  electron  - 
neutral  atmosphere  collisions. 

3.  Particle-in-Cell  (PIC)  beam  launch  and  propagation  model:  The  TRISTAN 
electromagnetic  and  relativistic  PIC  code  [Bmeman  et  ah,  1993]  was  applied  to  beam  injection 
from  a  spacecraft  to  asses  spacecraft  charging,  beam  degradation,  and  questions  of  stability. 

4.  Khazanov  plasma  transport  model:  Based  on  the  Boltzman  kinetic  equation  for  the  electron 
distribution,  the  Khazanov  4D  and  5D  plasma  transport  models  were  extended  to  relativistic 
energies.  The  time  dependent  behavior  of  relativistic  beam  flux  and  its  secondaries  in  the 
atmosphere-ionosphere-magnetosphere  was  modeled. 

The  first  two  are  discussed  in  the  following  section  with  the  latter  two  covered  in  the  next  two  sections. 

2.  The  Paraxial  Envelope  Equations  and  RelBAI  Model 

In  order  to  guarantee  a  successful  mission  in  which  the  heritage  of  an  accelerator  such  as  the  LINAC 
system  may  be  established,  it  is  useful  to  quantify  the  beam  propagation  characteristics  that  are  expected 
during  an  active  experiment  in  the  terrestrial  space  environment.  Specifically,  the  areas  of  interest 
covered  by  this  study  include: 

•  Analytical  Description  of  Beam  Propagation.  As  an  extension  to  the  study  by  Neubert  et  al 
[1996],  the  beam  envelope  equations  have  been  used  to  predict  beam  propagation  characteristics 
of  interest  (energy,  radius,  divergence)  as  a  function  of  altitude  over  a  variety  of  injection 
parameters,  including  beam  energy  (1  —  100  Me V),  divergence  (0.01  —  0.1  radians),  and 
relativistic  electron  beam  study  current  (0.1  -  10'‘  A).  Effects  of  beam-generated  electric  and 
magnetic  fields  were  included  in  the  analysis. 

•  Beam-Plasma  Interactions  (BPI).  When  high  current  beams  (I  >  100  A)  are  considered,  it  is 
appropriate  to  examine  the  interaction  of  the  beam  with  the  surrounding  space  plasma.  A  first- 
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order  analysis  of  the  BPI  was  performed  for  relativistic  electron  beams  (REBs)  propagating  in  a 
model  terrestrial  ionospheric  environment. 

•  Diagnostics  for  Active  Space  Experiments  with  REBs.  The  Relativistic  Beam-Atmosphere 
Interactions  (RelBAI)  model  has  been  constructed  in  order  to  provide  scientific  and 
technological  investigators  with  a  tool  to  quantify  physical  signatures,  such  as  optical  emissions 
and  ionization,  resulting  from  the  REB  precipitation  in  a  model  atmosphere.  This  model 
provides  a  means  of  predicting  the  intensity  of  these  products  that  would  aid  in  the  selection  of 
the  appropriate  instrumentation  for  detection  and  characterization  of  a  propagating  REB  during 
an  actual  test  in  space. 

This  section  presents  a  summary  of  the  findings  from  the  analysis  of  REB  dynamics,  stability,  and  the 
resulting  ionospheric  modification  during  propagation  in  a  model  terrestrial  upper  atmospheric  and 
ionospheric  environment.  References  to  published  technical  reports  are  noted  in  each  section  to  direct  the 
interested  reader  to  documents  detailing  this  research. 

2.1  Beam  Dynamics 

We  have  adopted  a  set  of  analytical  equations  that  describe  relativistic  electron  beam  propagation  in  the 
paraxial  approximation  (small  pitch  angle  and  narrow  divergence)  to  a  set  of  beams  that  are  propagating 
in  the  Earth's  near  space  environment.  The  model  has  application  to  naturally  occurring  space  science 
phenomena,  but  our  analysis  is  targeted  on  the  ideal  case  of  artificially  injecting  relativistic  elections 
beams  of  various  energies  (1-100  MeV),  divergence  values  (0.01  -  0.1  radians),  and  currents  (0.1  - 
lO"*  A).  These  beams  are  simulated  as  injected  downwards  along  a  constant  vertical  magnetic  field  (0.4 
G),  and  they  propagate  through  a  model  atmosphere  and  ionosphere  representative  of  a  night-time 
environment  with  moderate  solar  activity  (sunspot  number  ~  30).  The  model  atmosphere  and  ionosphere 
used  for  these  computations  is  detailed  in  Habash  Krause  [1997a].  The  radial  expansion  of  the  beam  as 
the  electrons  collide  with  the  atmospheric  neutrals  is  described  by  the  paraxial  ray  equation;  energy  loss 
due  to  elastic  and  inelastic  collisions  and  electron-neutral  bremsstrahlung  is  used  to  determine  pitch  angle 
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scattering  (emittance  growth)  of  the  beam.  The  complete  set  of  envelope  equations  are  given  as 
[Lawson,  19SS;  Humphries,  \99Q’,  Neubert  et  al.,  1996]: 
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where  Equations  2.1  through  2.6  respectively  describe  the  radial  expansion  of  the  beam,  the  emittance 
growth  due  to  all  processes,  the  emittance  growth  specifically  due  to  elastic  collisional  processes,  the 
energy  lost  by  the  beam  due  to  electron-neutral  collisions,  the  bremsstrahlung  energy  radiated  away  from 
the  beam,  and  the  beam  perveance. 


Sample  computations  were  performed  for  five  beam  energies  ranging  from  1.0  MeV  to  100  MeV.  All 
beams  carried  a  current  of  0.1  A  and  were  injected  with  an  initial  divergence  of  0.01  radians.  Figures 
2.1  through  2.3  respectively  show  the  radius,  mean-squared  divergence,  and  fractional  energy  loss, 
which  are  computed  at  each  altitude  from  injection  (140  km)  down  to  the  altitude  at  which  the  bulk  of  the 
energy  is  lost.  These  results  show  that  all  beams  maintain  a  slowly-varying  radius  until  they  encounter  a 
region  in  the  atmosphere  that  is  dense  enough  to  cause  significant  emittance  growth  via  electron-neutral 
collisions.  The  altitude  at  which  this  occurs  depends  on  the  initial  energy  of  the  beam,  which  is 
consistent  with  the  fact  that  the  interaction  cross  section  decreases  with  increasing  primary  beam  energy. 
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The  atmospheric  range  increases  monotonically  with  beam  primary  energy  {R(E)  =  3.247£’'  R  in 
kg/m-,  £  in  MeV),  but  it  is  found  to  be  independent  of  divergence  and  current.  The  range  values  agree 
well  with  those  computed  in  the  continuous  loss  approximation  [Berger  and  Seltzer,  1964],  but  they 
diverge  slightly  as  lower  energies  (£  <  5  MeV)  are  approached  (Figure  2.4).  See  Habash  Krause 
[1997a]  for  further  details  on  model  construction  and  computations.  Beam  propagation  following 
upward  injection  has  been  investigated  with  this  model  and  is  described  in  Habash  Krause  et  al  [2000]. 


Figure  2.1:  Beam  radius  as  a  function  of  altitude  for  downward  propagation.  Injection  current  and 
divergence  values  are  0.1  A  and  0.01  radians,  respectively.  Five  beams  are  examined  with 
injection  energies  E  -  [1,  5,  10,  50,  100]  MeV.  Inset  figure  contains  an  enlarged  view  of 
the  expansion  of  the  1,  5,  and  10  MeV  beams  as  they  enter  the  collisional  regime. 
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Figure  2.2:  Beam  divergence  as  a  function  of  altitude  for  downward  propagation.  Injection  current 
and  divergence  values  are  G.l  A  and  0.01  radians,  respectively.  Five  beams  are  examined 
with  injection  energies  E  -  [1,5,  10,  50,  100]  MeV.  Paraxial  limit  (0.78  radians)  is 
shown  for  reference. 


Figure  2.3:  Beam  fractional  energy  loss  per  electron  as  a  function  of  altitude  for  downward 
propagation.  Injection  current  and  divergence  values  are  0.1  A  and  0.01  radians, 
respectively.  Five  beams  are  examined  with  injection  energies  E  -  [1,  5,  10,  50,  100] 

MeV. 
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Figure  2.4:  Electron  beam  atmospheric  collisional  range:  Comparison  with  previously  published 
results. 


2.2  Beam  Stability 

Growth  rates  for  several  instabilities,  including  the  two-stream,  resistive  hose,  and  ion  hose,  were 
assessed.  The  equations  governing  the  time  rate  of  growth  of  the  instabilities  are: 
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where  equations  2.7  through  2.9  respectively  describe  the  characteristic  growth  periods  of  the  two- 
stream,  ion  hose,  and  resistive  hose  instabilities.  We  seek  these  values  because  if  the  growth  period  of 

an  instability  is  comparable  to  or  shorter  than  the  length  of  a  beam  pulse,  the  pulse  may  be  affected  by  the 
instability. 
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Additionally,  beam  susceptibility  to  Weibel  filamentation  was  examined  using  the  filamentation  stability 
factor,  given  as: 


2  2 
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where  ^  >  1  signifies  that  the  beam  is  stable  against  filamentation. 
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Using  the  same  model  atmosphere/ionosphere  as  those  employed  in  the  computation  of  beam  dynamics, 
instability  growth  periods  and  filamentation  stability  factors  were  computed  for  downward-injected 
beams  of  several  initial  energies  (1-10  MeV),  divergence  values  (0.01  -  0.1  radians),  and  currents 
(1  -  lO"*  A),  and  at  noon  and  midnight  conditions.  An  injection  altitude  of  350  km  was  selected  to 
investigate  the  worst-case  scenario  of  propagation  through  the  F  region.  It  was  found  that  within  this 
range  of  parameters,  all  of  these  beams,  when  pulsed  at  p,s  intervals,  are  stable  against  filamentation  and 
ion  hose,  but  low  energy,  high  cunent  beams  may  be  subject  to  the  resistive  hose  and  two-stream 
instabilities.  Sample  computations  appear  in  Figures  2,5  through  ZIO,  with  initial  beam  parameters 
given  in  the  captions.  The  ramifications  of  the  onset  of  any  of  these  instabilities  include  enhanced 
emittance  growth,  heating  and  ionization  of  the  ambient  plasma  environment,  and,  in  some  of  the  more 
esoteric  events,  current  multiplication.  These  factors  may  play  a  significant  role  in  the  description  of  the 
beam  dynamics.  If  so,  a  unified  approach  to  solving  the  dynamics  equations  with  these  effects  included 
would  be  necessary  for  a  complete  and  accurate  description  of  relativistic  electron  beams  propagating  in 
the  Earth’s  atmosphere.  A  detailed  treatment  of  the  instability  analysis  appears  in  Habash  Krause 
(1997b]. 


Figure  2.6:  Two-stream  instability  gro 


Figure  2.7:  Two-stream  instability  growth  period  for  nighttime  propagation;  E  =  10  MeV  and  <0‘>o 
=  O.I  rad. 


Figure  2.8:  Ion  hose  instability  growth  period  for  nighttime  and  daytime  propagation  of  all  space- 
charged  neutralized  beams. 
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2.3  The  RelBAI  Model 


The  RelBAI  model  was  constructed  for  the  purpose  of  providing  a  tool  by  which  investigators  interested 
in  detailed  specifications  of  products  of  the  relativistic  BAI  can  compute  quantities  of  interest  with 
minimal  computer  resources.  The  model  is  comprised  of  several  established  analytical  algorithms,  each 
of  which  is  designed  with  its  own  set  of  input  and  output  parameters.  The  idea  is  to  have  a  “black  box” 
model  that  allows  users  to  enter  beam  parameters,  environmental  models,  and  cross  section  Hata- 
quantities  of  interest  are  output  in  steady  state.  The  model  is  modular  in  nature,  and  input/output 
parameters  of  the  modules  are  wired  together  which  form  the  connections  of  the  modules  (Figure  2.11). 
This  arrangement  facilitates  specialization  of  the  code:  investigators  may  build  upon  the  model  by 
connecting  new  modules  at  the  appropriate  places.  As  a  demonstration  of  the  capabilities  of  this  model, 
several  quantities  of  interest  for  remote  sensing  purposes  were  computed,  including  total  ionization, 
bremsstrahlung  production  and  production  of  several  excited  states  of  neutral  and  ionized  molecular 
nitrogen  resulting  from  the  propagation  of  a  5  MeV  beam  through  the  atmosphere.  For  example, 
ionization  of  major  neutral  constituents  appears  in  Figure  2.12.  We  note  that  the  model  requires 
validation  with  experimental  data  once  the  SLINAC  is  operational,  but  under  this  contract,  we  were  at 
least  able  to  demonstrate  consistency  with  certain  acceptable  values  associated  with  electron  precipitation 
in  the  terrestrial  atmosphere.  For  example’  it  is  well-known  that  the  average  amount  of  energy  deposited 
per  ionization  event  in  the  atmosphere  is  ~35  eV,  and  our  results  are  comparable  to  this  number  (Figure 
2.13).  See  Hahash  Krause  [1998a]  for  a  detailed  treatment  on  the  construction  and  use  of  the  RelBAI 
model. 

In  the  next  section,  the  bremsstrahlung  fluxes  computed  from  this  model  are  used  to  predict  photon 
energy-dependent  spectra  incident  on  detectors  on  balloons  or  satellites  in  low  Earth  orbit. 
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RelBAI  Model  Block  Diagram 


Figure  2. 1 1 :  Block  diagram  of  the  RelBAI  model. 
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ical  electron  fluxes 


Ionization  by  fast  electrons,  E  =  5  MeV,  I  =  1 .0  A 


Figure  2.12:  Ionization  of  major  atmospheric  neutral  constituents  from  cascading  electrons  for  initial 
beam  energy  of  5  MeV. 

Deposited  Energy/Ionization  Event,  E  s  5  MeV,  I  =  1.0  A 


Figure  2.13:  Energy  deposition  per  ionization  event  from  cascading  electrons  for  initial  beam  energy  of 
5  MeV. 
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2.4  Remote  Sensing  of  Beam-Induced  Bremsstrablung 

Bremsstrahlung  flux  spectra,  differential  in  photon  energy  and  which  are  incident  on  detectors  at  ball(X)n 
and  Low-Earth  Orbit  (LEO)  satellite  altitudes,  are  computed.  Attenuation  of  radiation  by  the  atmosphere 
is  included,  and  a  mass  attenuation  coefficient,  dependent  on  photon  energy,  incorporates  intensity  loss 
mechanisms  of  interest  for  the  photon  energy  regime  under  consideration.  These  scattering  mechanisms 
include  photoelectric  effect,  Rayleigh  scattering,  Compton  scattering,  electron-positron  pair  production, 
and  photonuclear  interaction.  The  geometry  for  computations  of  the  bremsstrahlung  flux  measured  by 
imagers  aboard  satellites  in  LEO  is  shown  in  Figure  2.14.  The  attenuated  bremsstrahlung  measured  by 
satellite  at  large  horizontal  distances  from  the  beam  column  is  larger  by  several  orders  of  magnitude  than 
that  measured  by  the  satellite  detector  at  small  horizontal  displacements  (Figure  2.15).  A  balloon  detector 
may  intercept  more  bremsstrahlung,  by  several  orders  of  magnitude,  as  in  the  situation  examined  here, 
than  the  downward-looking  satellite  detector  if  the  balloon  detector  is  oriented  with  its  look  direction 
coincident  with  the  beam  source  line  of  sight  (compare  Figures  2.15  and  2.16).  However,  it  may  be 
impractical  to  have  sufficient  knowledge  of  the  location  of  the  bremsstrahlung  source  before  the  balloon 
detector  can  be  reoriented  to  look  directly  at  the  beam.  Perhaps  a  balloon  payload  with  x-ray  sensors 
looking  in  each  of  the  four  cardinal  directions  would  be  necessary.  See  Habash  Krause  [1998b]  for 
details  on  these  computations. 
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Integrated  bremsstrahlung  flux  at  balloon  altitudes,  with  attenuation 


Figure  2.16:  Integrated  bremsstrahlung  flux  incident  on  a  detector  at  balloon  altitudes.  Atmospheric 
attenuation  is  included 

2.5  RelBAI  Summary  and  Continuing  Work 

The  research  completed  under  the  contract  has  laid  the  foundation  for  continuation  of  theoretical  studies 
necessary  to  plan  for  the  heritage  launch  of  the  SLINAC  5  MeV  electron  beam.  The  next  major  step  in 
the  development  of  this  research  was  accomplished  on-site  at  AFRI7VSBS  with  the  completion  of  a 
relativistic  electron  auroral  model.  This  model  contains  the  time-dependent  chemistry  necessary  to 
compute  optical  emissions  as  a  function  of  altitude.  It  was  shown  that,  for  5  MeV  electron  beams 
injected  directly  downward  into  the  Earth’s  atmosphere  along  a  constant,  vertical  magnetic  field,  the 
427.8  nm,  391.4  nm,  and  337.1  nm  molecular  nitrogen  lines  are  strong  in  intensity,  with  production 


19 


rates  several  orders  of  magnitude  greater  in  than  those  of  the  O  lines  examined  (557.7  nm  and  630.0- 
636.4  nm  doublet).  It  was  shown  that  the  production  of  337.1  nm  is  insensitive  to  compositional  change 
and  has  a  quenching  height  lower  in  altitude  than  the  propagation  depth  of  a  5  MeV  electron  beam,  and 
thus  the  signature  may  be  suitable  as  an  indicator  of  electron  flux  for  beams  of  comparable  energy.  The 
ratio  of  427.8  nm  to  391.4  nm  emissions  was  found  to  be  relatively  insensitive  to  compositional 
changes,  and  the  ratio  varies  with  altitude  at  lower  altitudes,  suggesting  that  it  may  suitable  for  inference 
of  characteristic  beam  energy  for  MeV  electron  beams  (see  Habash  Krause  et  al.  [1999]  for  details  on  the 
relativistic  electron  auroral  model). 

3.  3-D  Relativistic  Particle-in-Cell  (PIC)  Simulations  of  The  Relativistic 
Beam  Injection  Process 

3.1  Background 

PIC  simulations  of  relativistic  electron  beam  dynamics  have  been  undertaken  to  investigate  the  beam 
injection  process  from  a  spacecraft  in  order  to  give  realistic  initial  conditions  of  beam  parameters  for 
large-scale  models  of  beam  propagation  in  the  Earth’s  atmosphere-ionosphere-magnetosphere  [Neubert  et 
al.,  1996;  Kazanov  et  at.,  1999].  In  particular,  the  following  questions  have  been  addressed: 

•  The  proper  handling  of  spacecraft  charging  effects  on  the  beam  in  a  3D  code.  This  is  an  important 
issue  because  beams  will  often  not  escape  a  spacecraft  environment  with  the  full  energy  because  of 
electric  charging  of  the  spacecraft. 

•  Kinetic  effects  influence  beam  propagation  in  the  plasma.  PIC  simulations  are  done  to  help 
complement  the  analytical  beam  dynamics  analysis,  in  particular  to  address  the  question  of  beam 
neutralization  and  beam  stability. 

The  beam  pulses  under  investigation  are  those  of  the  linac  being  tested  for  flight  by  the  Air  Force  [Jorr, 
1993].  The  parameters  are  shown  in  Table  3.1. 
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Table  3 . 1 :  Relativistic  beam  pulse  parameters 


Ib 

Eb 

rb 

Td 

ao 

e 

nb° 

(opb° 

O.I  A 

5  MeV 

4  ps 

0.001 

1.5  mm 

0.018  cm-rad 

EESQIjQIH 

EEniiysiinii^ 

4  is  the  beam  current,  the  beam  energy,  the  beam  pulse  length,  the  duty  cycle,  a„  the  beam 
diameter  at  injection,  e  the  beam  emittance,  the  density  at  injection,  and  the  beam  plasma 
frequency  at  injection:  where  y  is  the  relativistic  factor. 

Simulations  were  first  performed  in  2D,  but  it  was  soon  realized,  that  3D  should  be  attempted  within  the 
framework  of  the  limited  effort  available  in  this  research.  The  3D  simulations  are  performed  using  the 
TRISTAN  electromagnetic  and  relativistic  particle  code  [Buneman  et  al.,  1993].  The  code  uses  local 
updates  of  the  fields  from  particle  motions,  rather  than  Poisson’s  equation.  While  this  scheme  makes  the 
code  fast,  it  is  required,  that  the  experimental  conditions  are  described  self-consistently  within  the 
simulation  domain.  It  is  important,  therefore,  to  incorporate  the  beam  source  and  the  spacecraft  in  the 
simulation  domain  with  care. 

Code  developments  during  this  study  include: 

•  Self-consistent  spacecraft  charging  due  to  beam  pulse  emission  has  been  included,  and  the 
appropriate  scaling  laws  between  simulations  and  the  real  world  for  spacecraft  charging  studies  have 
been  identified. 

•  The  option  of  changing  to  a  moving  frame  which  follows  the  beam  pulse  has  been  included,  but  is  not 
used  here. 

Conclusions  obtained  during  the  study  can  be  summarized  as: 

•  The  spacecraft  charging  levels  generated  during  single  pulse  injections  are  small  compared  to  the 
beam  energy  and  a  pulse  will  escape  the  spacecraft  loosing  only  a  small  fraction  of  its  energy  -  even 
in  vacuum.  During  repeated  pulsing,  the  effective  current  is  0.1  mA,  which  is  small  compared  to  the 
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thermal  electron  current  to  a  spacecraft  in  the  ionosphere.  Charging  rates  are  expected  to  be  the  same 
as  for  non-relativistic  beams,  i.e.,  to  follow  the  formula  of  Parker  and  Murphy  [1967],  with  an 
enhancement  factor  of  2.5  for  low  earth  orbit  spacecraft  as  demonstrated  in  the  TSS-IR  mission 
[Thompson  et  al.,  1998].  Significant  charging  (relative  to  the  beam  energy)  is  expected  only  when  the 
ambient  electron  density  is  below  lO^m’^. 


•  Relativistic  beams  seem  to  propagate  with  higher  stability  than  keV  beams  studied  in  the  past, 
probably  because  of  the  larger  beam  electron  velocity,  which  reduces  space  charge  effects,  and  larger 
relativistic  electron  mass.  Even  over-dense  beams,  where  the  beam  density  at  injection  exceeds  the 
plasma  density  by  several  order  of  magnitudes,  propagate  well. 

•  A  new  phenomenon  has  been  identified  for  the  relativistic  regime,  which  is  “ion-focusing”,  that  may 
be  of  importance  for  beam  currents  over  10  A,  but  will  not  affect  the  linac  beam  at  lower  currents. 
In  the  ion-focused  regime,  a  beam  is  partially  space  charge  neutralized.  The  beam  self-forces  tend  to 
focus  the  beam,  enhancing  the  process  of  pushing  ambient  electrons  out  of  the  beam  region,  which 
leads  to  further  beam  neutralization. 

3.2  Spacecraft  Charging 

3.2.1  The  Plasma  Response  Time 

The  spacecraft  potential  will  initially  increase  with  time,  causing  beam  electrons  to  escape  the  spacecraft 
environment  with  decreasing  energy.  The  ambient  plasma  will  react  and  attempt  to  supply  an  electrically 
neutralizing  return  current.  In  the  process,  an  electric  potential  sheath  is  established  in  the  plasma  around 
the  spacecraft,  causing  ambient  electrons  incident  on  the  sheath  edge,  to  be  drawn  onto  the  spacecraft. 

In  steady-state,  the  sheath  edge  is  located  at  a  distance  from  the  spacecraft  which  allows  a  return  current 
that  equals  the  injected  beam  current.  For  ionospheric  sounding  rockets,  the  relationship  between  the 
beam  current  and  the  potential  has  been  found  to  follow  the  formula  of  Parker  and  Murphy  [1967], 


22 


which  include  effects  of  ambient  magnetic  fields  [Myers  et  al.,  1989].  Later  investigations  of  data  from 
the  NASA  TSS-IR  mission  showed,  that  the  functional  dependence  between  current  and  potential  for 
LEO  satellites  also  followed  Parker  and  Murphy,  but  with  an  enhancement  factor  of  about  2.5, 
presumably  caused  by  the  high  orbital  motion  of  LEO  spacecraft  [Thompson  et  al.,  1998].  The  good 
agreement  between  theory  and  experiments  was  obtained  for  beams  with  energies  in  the  few  keV-range. 
The  relativistic  beams  investigated  here  are  not  expected  to  behave  any  differently,  and  we  conclude 
based  on  the  work  by  Myers  et  a/.  [1989],  that  significant  charging  during  repeated  pulsing  of  the  linac 
beam  will  begin  only  when  the  ambient  density  drops  below  10’m‘\ 

It  is  worthwhile  pointing  out,  the  beam  experiments  in  space  and  related  beam  simulations  studies  of  the 
past  considered  beam  pulses  with  pulse  lengths  of  the  order  of  milliseconds  or  more  [Neubert  and 
Banks,  1992].  Studying  the  linac  p,s-duration  pulses  corresponds  to  investigating  the  beam  front  in  these 
experiments.  The  key  parameter  for  the  charging  state  during  pulse  emissions  is  the  plasma  response  time 
which  is  of  the  order  of  l/co^.  The  ambient  plasma  density  for  which  is  approximately  n^  =  l  y, 

10’ m‘^.  For  comparison,  f-peak  ionospheric  densities  typically  range  from  10'“  (night)  to  10'^  m'^  (day). 
In  the  ionosphere,  steady  state  conditions  are  then  expected  to  be  reached  during  a  linac  pulse  emission. 
Regions  of  low  plasma  density  are  found  both  below  and  above  the  ionosphere.  Below  the  ionosphere, 
the  beam-spacecraft  system  is  dominated  by  collisional  interactions  with  the  atmosphere,  which  lead  to 
the  creation  of  a  substantial  flux  of  secondary  electrons  [Neubert  et  al.,  1990],  a  regime  not  considered 
here.  However,  well  above  the  ionosphere,  in  the  magnetosphere,  the  plasma  density  may  be  so  low  as 
to  create  conditions  for  the  transient  state. 

From  the  above  discussion,  it  is  clear,  that  the  first  parameter  to  scale  properly  in  the  simulation  is  the 
ratio  of  the  beam  pulse  time  to  the  plasma  response  time,  =  t/t^. 

3.2.2  The  spacecraft  charging  rate 

A  second  process  to  scale  properly  is  the  charging  experienced  by  a  spacecraft  during  a  pulse  injection. 
The  charging  rate,  if  large,  may  severely  affect  the  beam  pulse  escape  energy  from  the  spacecraft  sheath 
and  modify  the  beam  dynamics.  The  charging  rate  of  a  capacitor  is: 


23 


^-L 

dt~  C 


(3.1) 


where  Vis  the  voltage,  / the  current,  and  C the  capacitance.  The  capacitance  of  a  conducting  sphere  with 
radius  r  in  vacuum  is: 

C  =  4nej  (3  2) 

With  r  =  1  m  and  /  =  I  A,  the  charging  rate  is  9  kV  p,s'‘. 

Assuming  that  a  1-m  radius  sphere  is  a  reasonable  approximation  of  a  spacecraft  for  charging 
considerations,  it  follows  that  the  beam  system  simulated  here,  with  4  =  0.1  A  and  =  4  ^is,  under  the 
transient  condition  will  charge  to  a  maximum  of  3.6  kV.  The  potential  energy  of  electrons  on  the 
spacecraft  at  the  end  of  a  pulse  injection  is  then,  even  in  worst  case,  small  relative  to  the  beam  accelerator 
energy,  but  large  relative  to  the  ambient  plasma  thennal  energy.  In  the  simulations,  the  scaling  must  be 
such  that  the  relative  magnitude  of  these  entities  are  maintained. 


3.2.3  Scaling  of  spacecraft  charging 

To  properly  compare  the  simulation  world  and  the  real  world,  one  needs  to  compare  the  rate  of  change  of 
the  nonmlized  potential 


(3.3) 


With  (3. 1-3.2),  (3.3)  becomes 


d4>  _  ql 
dt  Anejnc^r 


The  important  relationship  to  establish  in  the  simulation  is  the  number  of  beam  electrons  to  inject  in  one 
time  step,  such  that  the  electric  charging  potential  energy  experienced  during  one  beam  pulse  injection  is 
small  compared  to  the  accelerator  energy,  but  large  compared  to  the  ambient  plasma  energy.  The  number 

depends  on  simulation  values  of  the  electron  mass  and  charge  and  the  number  density  of  the  ambient 
plasma. 
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The  beam  current  at  the  accelerator  can  be  given  in  terms  of  the  number  of  electrons  ejected  per  unit  time 
Na': 


h  =  qK 

With  (3.4-3. 5),  the  charging  state  can  now  be  described  by: 

AO  =  ^  ”  .  At. 

Ansjnrc 


(3.5) 


(3.6) 


or: 

AO  = - ,  Af 

ATtn^rCX  p 


(3.7) 


For  A/  =  ar^,  where  a  is  the  fraction  of  a  beam  pulse  injected,  we  can  write  the  charging  state  as: 


AO  = 


ATOt^rc^ 


(3.8) 


Equation  3.8  is  used  for  deciding  on  the  simulation  parameters.  The  simulation  domain  is  on  a  3D  grid 
with  grid-spacing  A=l,  and  time  step  At  =  1.  The  velocity  of  light,  c  =  0.5,  such  that  a  relativistic 
electron  will  move  approximately  0.5A  in  one  time  step.  The  first  simulation  grid  used  is  256A  x  128A  x 
128A,  with  the  beam  pulse  propagating  in  the  positive  x-direction. 


The  spacecraft  is  as  a  cube  of  size  lOA  x  lOA  x  lOA,  centered  in  y  and  z  and  at  a:  =  64A,  such  that  the 
approximate  equivalent  spherical  spacecraft  radius  is  r  =  5.  For  simulations  of  a  complete  beam  pulse  (a 
=  1),  the  beam  pulse  duration  is  chosen  to  be  250  At,  which  allows  both  the  head  and  tail  of  the  pulse  to 
be  resolved  properly.  The  ambient  simulation  density  is  n  =  1  for  the  steady  state  conditions,  which  is  a 
low  number  to  reduce  computational  effort,  but  not  to  low  to  resolve  kinetic  plasma  effects.  The  transient 
state,  where  no  significant  return  currents  are  generated  during  the  time  of  the  pulse  injection,  will  be 
simulated  using  injection  into  a  vacuum. 


For  steady-state,  we  chose  a  plasma  with  approximately  five  ambient  plasma  periods  during  one  beam 
pulse  injection,  or  more  specifically,  xJT^  =  5.17,  giving  x^^  =  32.5.  In  the  simulation,  ions  are 
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assumed  stationaiy.  Wid,  the  code  employed,  one  need  not  include  ions  at  all  in  this  case,  thus  further 
reducing  the  computational  load. 


With  these  numbers,  (3.8)  becomes: 


A<P  =  a — ~ 
3.72 


(3.9) 


When  deciding  on  the  beam  panicle  injection  rate,  the  requirement  is  to  maintain  the  relative  oniering  of 
normalized  plasma  energy,  beam  energy  and  final  spacecraft  electric  potential  energy  for  vacuum.  From 
(3.9)  we  conclude  that  10  electrons  pr.  dme  step  seems  a  reasonable  injection  rate.  With  a  beam  cross 
section  at  injecdon  in  the  simulations  of  4A=  and  a  simuladon  velocity  of  0.5,  the  beam  density  at 
injection  is  then  5,  or  5  times  the  ambient  density  (ii,7«,=5).  The  normalized  energies  for  /V,  ’=  10  is 
shown  in  Table  3.2.  With  our  choice  of  parameters,  the  relative  order  is  maintained  in  the  simulation. 
However,  the  escape  energy  in  vacuum  will  be  reduced  to  7. 1,  which  is  more  than  in  the  real  world.  As 
will  be  clear  from  the  simulations,  this  is  not  a  problem  for  the  physics  under  investigation. 


Table  3.'2:  normalized  energies  representative  for  an 
experiment  in  the  ionosphere  and  in  the  simulation 


-  at - -  — 

real  worid 

simulation 

plasma  enerev 

IF 

2xTF  " 

beam  energy 

9.8 

9.8 

S/C  potential 

0.007 

^.7 

3.2.4  Test  of  Spacecraft  Charging  Rates 

The  scaling  law  and  the  implementation  of  a  spacecraft  in  the  simulation  has  been  tested  by  injection  of  a 
beam  into  a  vacuum  with  the  parameters  given  In  the  past  subsection.  The  normalized  potential  of  the 
spacecraft  is  determined  each  time  step  by  integrating  the  radial  electric  field  component  along  y  from  the 
j=0,bounda.y  and  to  the  center  of  the  spacecraft.  Particles  hitting  the  spacecraft  are  arrested  and  written 
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off,  and  electric  fields  from  surplus  charge  are  conducted  freely  within  the  spacecraft  structure.  Thus, 
charge  from  the  accelerator  is  redistributed  with  the  speed  of  light  within  the  structure. 

Figure  3.1  shows  the  spacecraft  potential  as  a  function  of  time  (full  line)  and  the  corresponding  analytical 
prediction  for  vacuum  given  by  (3.4)  (dashed  line).  On  the  top  is  shown  the  vacuum  case  and  on  the 
bottom  the  plasma  case.  The  electric  field  forming  around  the  spacecraft  and  the  injected  beam  is  being 
established  with  the  speed  of  light.  The  charging  of  the  spacecraft  is  thus  lagging  the  analytical  estimate. 
However,  the  charging  rates  for  vacuum  shown  on  the  left  compare  veiy  well,  and  the  spacecraft  retains 
it  potential  after  termination  of  pulse  injection  (not  shown).  For  injection  into  a  plasma  {n°ln=5),  the 
spacecraft  is  clearly  well  neutralized  by  the  plasma  current.  The  time  scale  is  sown  in  units  of  the  plasma 
period  for  the  plasma  case  on  the  bottom  panel. 

Thus,  we  are  satisfied,  that  the  incorporation  of  a  spacecraft  and  accelerator  into  the  TRISTAN  code  is 

done  properly,  and  that  the  correct  scaling  laws  between  the  real  world  and  the  simulation  world  have 
been  identified. 

3.3  Beam  Dynamics 

3.3.1  A  First  Look  at  Beam  Dynamics 

A  first  cut  at  simulations  of  a  complete  beam  pulse  is  done  with  the  parameters  given  in  the  previous 
section  for  the  transient  condition  with  injection  into  a  vacuum  and  for  the  steady  state  conditions  with 
^bp  =  5.17  in  "hi  =5).  Results  are  given  at  a  point  just  after  the  beam  is  turned  off  at  t  =  256At,  where 
the  beam  pulse  has  been  injected  in  the  interval  t  =  0-250At.  Figure  3.2  shows  the  position  and  the 
energy  of  the  beam  electrons  for  the  vacuum  case  (left)  and  the  ambient  plasma  case  (right).  The  beam 
has  been  injected  at  45“  pitch  angle  with  8“  angular  spread.  The  non-relativistic  gyro  frequency  equals  the 
non-relativistic  plasma  frequency  for  the  plasma  case. 
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Figure  3.1: 


0  1  2  3  4  5 

Time  [Tpe] 


Spacecraft  potential  as  a  function  of  time.  Top:  beam  into  vacuum;  Bottom:  beam  into 
plasma  (nb°lnp=5).  Full  line:  simulation,  dashed  line:  vacuum  theory  (equation  3.4). 
Time  is  measured  in  plasma  periods  for  the  plasma  case. 
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Figure  3.2:  Beam  position  and  energy  at  t=256At  {tlTpe=52  for  the  plasma  case).  Beam  injection 
pitch  angle  6  =  45°  with  angular  spread  40  =  8°;  beam  density  at  injection.  Left:  beam 
into  vacuum;  right:  beam  into  plasma  («/  trip  =  5).  An  ambient  magnetic  field  is  directed 
along  the  x-axis. 
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Looking  first  at  the  beam  position,  it  is  noted,  that  the  beam  in  vacuum  has  undergone  significant  radial 
expansion,  which  primarily  is  due  to  the  initial  spread  in  the  beam  injection  angle,  while  in  the  plasma 
case  the  beam  is  much  more  focused.  From  the  plots  of  beam  energy  as  a  function  of  position,  it  is  clear, 
that  in  the  vacuum  case,  the  beam  experiences  a  gradual  decrease  in  energy,  which  corresponds  to  the 
loss  of  escape  energy  because  of  spacecraft  charging  consistent  with  the  charging  shown  in  Figure  3.1 
(top),  for  the  same  simulation  parameters.  The  energy  loss  is  smooth  and  the  beam  pulse  undergoes  no 
unstable  oscillations.  In  the  plasma  case  it  is  clear  that  the  beam  pulse  is  subject  to  strong  oscillations, 
leading  both  to  higher  and  much  lower  energies  than  the  injection  energy.  The  oscillations  grow  with 
time  in  the  plasma,  being  more  severe  in  the  front  of  the  pulse  than  in  the  rear.  It  is  also  seen  from  the  tail 
of  the  pulse,  that  the  escape  energy  in  the  plasma  is  close  to  the  accelerator  energy,  i.e.,  the  spacecraft  is 
well  neutralized. 


Because  the  beam  density  relative  to  the  ambient  density  clearly  is  an  important  parameter,  we  must  look 
more  closely  at  the  injection  process,  and  carefully  simulate  the  beam  radial  expansion  process.  But  first 
we  take  a  look  at  an  useful  analyticaf  description  of  the  beam  expansion  problem. 

3.3.2  The  Paraxial  Envelope  Equation 

The  radial  expansion  of  the  beam  is  described  by  the  paraxial  envelope  equations  [Lawson,  1988, 
Humphries,  1990].  Considering  only  terms  relevant  for  the  present  problem  of  a  relativistic  electron 
beam  in  a  magnetic  field,  the  equations  reduce  to: 


3  o2  2 

a  a  pc 
(1)  (2)  (3) 


(3.10) 


where  a”  is  the  second  derivative  of  the  beam  radius  with  respect  to  the  distance  along  the  beam  axis, 
and  is  the  beam  electron  Larmor  frequency  (=-qB/2ymJ.  K  is  the  generalized  perveance: 


2/ 


(3.11) 


where  the  normalized  current  and: 
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(3.12) 


(=17000  A) 


The  perveance  is  the  sum  of  beam-generated  forces  from  the  non-neutralized  space  charge  within  the 
beam,  that  attempts  to  expand  the  beam,  and  the  beam  self  magnetic  field,  which  tends  to  contract  the 
beam.  The  parameter/is  the  degree  of  beam  charge  neutralization  by  the  ambient  plasma,  and  it  takes  on 
values  between  0  and  1.  A  completely  charge  neutralized  beam  will  have  f=l,  in  which  case  K  is 
negative.  It  is  interesting  to  look  at  the  perveance  values  in  the  limits  ^0,1  for  non-relativistic  (}<=1)  and 
relativistic  (^1)  beams.  At  these  limits  we  have: 


-21  <K<  21/0^  non-relativistic 

(3.13) 

-21Iy  <K <2 I/y^  relativistic 

(3.14) 

For  non-relativistic  beams,  K  is  likely  to  be  positive,  while  for  relativistic  beams,  K  will  often  be 
negative.  From  (3.10),  one  can  see  that  the  sign  of  K  is  important  for  the  propagation  characteristics. 
Even  with  no  ambient  magnetic  field,  relativistic  beams  with  K  negative  have  the  possibility  of 
propagation  with  the  radius  maintained  below  a  maximum  value.  In  this  case,  the  positive  emittance  tema 
balances  the  negative  perveance  term.  Moving  from  ndn-relativistic  beams  to  relativistic  ones,  one  then 
has  to  consider  the  new  phenomena  of  beam  propagation  in  the  ion-focused  regime. 


3.3.3  Simulation  of  Radial  Expansion 

The  radial  expansion  of  the  beam  is  very  important  to  understand  because  the  stability  of  the  beam 
depends  on  the  ratio  of  the  beam  and  the  ambient  plasma  density.  For  low  beam  density,  the  interaction  is 
weak.  To  better  understand  the  relative  importance  of  the  three  terms  in  the  envelope  equation,  and  to 
help  assure  that  they  are  scaled  properly  to  the  simulations,  values  of  the  terms  at  injection  are  computed 
for  a  range  of  beam  currents  and  shown  in  Table  3.3,  The  table  shows  values  of  the  perveance  term  for/ 
=  0  {K,^  and  for/=  I  (AT,).  First  it  is  noted,  that  at  injection,  the  emittance  term  (2)  dominates.  It  can  be 
shown,  that  initially,  the  beam  expands  with  the  emittance  angle  Q^s/a^.  In  case  of  a  charge  neutralized 
beam,  the  next  term  to  dominate  when  the  radius  increases  is  the  perveance  term.  Finally,  if  the  beam 
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continues  to  expand,  the  ambient  magnetic  field  will  eventually  bring  the  radius  back  in  again.  The  last 
columns  show  the  radii  where  the  magnitude  of  the  terms  equal  each  other. 


The  current  values  shown  in  the  table  range  from  the  beam  current  of  0.1  A  and  to  100  A.  As  will  be 
shown  later,  it  is  only  for  the  highest  values,  that  ion-focusing  may  become  an  issue,  at  least  predicted  by 
the  envelope  equations.  With  the  range  shown,  we  thus  explore  the  regime  from  below  to  above  tfie 
currents,  where  this  may  become  an  issue  in  a  space  experiment,  and  compare  these  analytical  predictions 
with  simulations,  where  plasma  kinetic  effects  are  accounted  for.  An  important  issue  here  is  to  what 
extent  a  beam  is  space  charge  neutralized.  This  is  what  a  particle  code  can  elucidate. 


Magnitude  of  terms  in  Equation  3.10  at  injection  in  ionosphere  and  in  simulation. 
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3.4 

To  simulate  the  radial  expansion  of  the  beam,  injection  is  at  0®  pitch  angle  along  x.  This  allows  a  smaller 
dimension  perpendicular  to  the  beam  and  the  ambient  magnetic  field,  and  the  simulation  domain  is  now 
1024A  X  64  A  X  64  A.  The  pulse  length  is  now  4  times  as  long  (1000  At)  and  the  electric  charge  one 
quarter  of  the  one  used  previously,  such  that  one  pulse  length  still  equals  5.17  Tpe  for  the  plasma  cases. 
The  gyro  radius  is  maintained  unchanged  to  limit  the  radial  excursions,  and  thus  the  gyro  frequency  is  4 
times  the  one  previously  used. 


The  terms  of  the  paraxial  envelope  equation  for  the  simulations  are  also  shown  in  Table  3.X  In  the 
simulation,  the  emittance  term  also  dominates,  except  for  the  highest  current  values  investigated.  The 
perveance  terms  scale  almost  directly  between  the  real  world  and  the  simulation  world.  The  emittance 
terms  are  generally  smaller  in  the  simulations  because  it  is  impossible  to  completely  resolve  a  small  beam 
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radius  and  at  the  same  time  contain  a  full  pulse  length  within  the  simulation  domain.  Recall,  that  the  beam 
radius  is  1.5  mm,  and  the  pulse  length  is  1.2  km.  In  the  simulation,  the  radius  is  therefore  relatively 
larger  than  in  the  real  world.  The  ambient  magnetic  field  term  (3)  is  also  the  smallest  in  the  simulation, 
except  for  the  lowest  current  case. 

Equation  3.10  has  been  solved  numerically  with  the  real  world  current  values  of  Table  3.3  and  results  are 
shown  in  Figure  3.3.  On  the  top  is  shown  normalized  radius  as  a  function  of  distance  along  the  beam  for 
the  various  current  cases  and  on  the  bottom  the  same  for  the  simulated  currents.  The  thick  line  represents 
the  solutions  for  all  current  cases  with  K=Ko,  and  for  the  lowest  current  cases  for  K=K,.  For  these 
situations  the  perveance  term  is  so  small  as  to  have  no  impact  on  the  solution.  Only  as  currents  increase, 
does  the  perveance  term  begin  to  influence  the  solution,  but  only  for  the  case  of  a  charge-neutral  beam 
(K=K,).  For  the  real  beam,  the  1.2  km  pulse  length  is  shown  in  a  typical  ionospheric  magnetic  field, 
which  corresponds  to  about  0.5  of  a  oscillation.  For  the  simulation,  the  magnetic  field  is  relatively 
higher,  and  about  3.5  oscillations  are  shown  for  the  cases  represented  by  the  thick  curve.  Figure  3.3  ’ 
demonstrates  that  the  simulations  also  cover  the  transition  region  from  low  currents,  where  the  perveance 
term  is  unimportant,  to  the  high  current  region,  where  ion-focusing  may  become  important. 

Simulation  results  are  shown  in  Figures  3-.4-3»6  for  beam  in  vacuum  (n/=:l),  and  beams  in  plasma 
=5,  50).  On  the  left  panels,  the  beam  electron  positions  (top)  and  energies  (bottom)  are  shown 
just  after  the  end  of  the  beam  pulse  emission  at  t/T^-53  and  on  the  right  panels  at  On  the  plots 

of  the  beam  electron  positions,  the  solutions  to  3.10  are  shown  for  K^Kq  (top  curves)  and  K-K, 
(bottom  curves)  (for  the  vacuum  case  both  curves  are  for  K-Kg  as  the  beam  must  be  completely  non- 
neutralized).  The  beam  envelope  modulation  seen  in  the  simulations  and  the  analytical  solutions  for  a 
non-neutralized  beam  are  in  very  good  agreement  in  the  first  2  figures,  and  mainly  reflects  the  gyration  of 
electrons  in  the  presence  of  the  ambient  magnetic  field.  In  the  last  case  shown  in  Figure  3.6,  however, 
the  picture  becomes  much  more  complex.  This  is  also  the  case,  where  a  negative  perveance  may  play  a 
significant  role.  This  term  pinches  portions  of  the  beam  front  enough  to  push  the  ambient  electrons  aside, 
and  allow  segments  of  the  beam  to  propagate  as  partially  charge  neutralized.  At  the  same  time,  the  beam 


is  undergoing  strong  beam-plasma  interactions  with  spatial  modulation  lengths  that  are  shorter  than 
predicted  by  3.10.  The  results  indicate,  that  while  the  analytical  considerations  point  to  the  beam  forces 
as  instrumental  in  establishing  conditions  for  beam  pinching  and  propagation  in  the  ion  focused  regime, 
beam-plasma  interactions  dominate  the  propagation  aspects  by  establishing  more  complicated  plasma 
structures  in  the  beam  region,  which  in  turn  determines  a  rather  complicated  functional  dependence  of  the 
parameter /  with  space  and  time. 


Figure  3.3:  Beam  envelopes  for  different  beam  currents.  Top:  real  beams,  4  =  1,  30,  60,  100  A.  Thick 
curve  is  for  non-neutralized  beams,  where  Ko  is  small  for  current  values  investigated,  as  well  as 
for  neutralized  beam  at  the  lowest  current,  where  K,  is  also  small.  Dashed  curves  are  for 
neutralized  beams,  where  increased  currents  have  increased  pinching  effect  on  beam.  Bottom: 
same  for  simulated  beams  with  ^‘=1,  10,  100.  For  real  beams  (top),  distance  along  beam  (x) 
shown  corresponds  to  length  of  a  4  ps-beam.  Envelope  of  thick  curve  represents  magnetic 
field  focusing  and  distance  shown  corresponds  to  app.  0.5  oscillation  length.  In  simulation,  a 
relatively  higher  magnetic  field  is  chosen  to  limit  beam  radius.  Here,  thick  curve  traces 
approximately  3.5  magnetic  field-induced  oscillations. 
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Figure  3.4:  Beam  in  vacuum  vv'ith  magnetic  field  directed  along  x.  Beam  density  is  «/,  -5,  magnetic  field 
X2/4/=5.6xl0'\  and  beam  pitch  angle  at  injection  0  =  0",  with  A0=8°.  Beam  pulse  length  t*- 
10004/.  Left:  /=  10244/;  Right:  /=  16244/.  Top:  beam  position  in  the  acy-plane.  Horizontal 
scale  is  highly  compressed  as  seen  in  distortion  of  spacecraft.  Full  curve  is  solution  to  beam 
envelope  Equation  3.10  for  non-neutral ized  (^0)  beam;  Bottom:  normalized  beam  energy  as 
a  function  of  x. 
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Figure  3.5:  Beam  into  plasma  {n^°lnp=5)  with  a  magnetic  field  directed  along  the  x-axis  (i2t4/=5.6xl0’^). 

Beam  pitch  angle  at  injection  d  =  0°,  with  Ad  =  8".  Beam  pulse  length  Ti=  1000^/.  Left: 
/=  10244/  =  5.3);  right:  /=16244/  =  8.4).  Top:  beam  position  in  x>'-plane; 

bottom:  beam  energy  as  a  function  of  x.  Full  curve  is  solution  to  beam  envelope  Equation 
3.10  for  non-neutralized  (f=0,  top  curve)  and  neutralized  (f=l,  bottom  curve)  beam. 
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Figure  3.6:  Beam  into  plasma  (n^Vn^^SO)  with  a  magnetic  field  directed  along  the  x-axis 
(i2^./d/=5.6xl0'^).  Beam  pitch  angle  at  injection  9  =  0°,  with  Ad  -  8°.  Beam  pulse  length 
r,=  1000^/.  Left:  f=10244/  (t/T^  =  5.3);  right:  ^=162421/  (//I^  =  8.4).  Top:  beam 
position  in  Ay-plane;  bottom:  beam  energy  as  a  function  of  x.  Full  curve  is  solution  to 
beam  envelope  Equation  (3.10)  for  non-neutralized  (f=0,  top  curve)  and  neutralized  (f-1, 

bottom  curve)  beam. 

To  expand  further  on  this  point.  Figure  3.7  shows  the  electron  density  along  the  x-axis  through  the 
center  of  the  beam  as  a  function  of  time.  On  the  left  plots  are  shown  the  case  of  n,7n,  =5  and  on  the  right 
n;/n^  =50.  Top  plots  are  the  beam  densities,  and  bottom  plots  the  ambient  densities.  The  plasma  is 
pushed  out  of  the  beam  region  at  the  nodes,  where  the  electrons  after  one  gyration  concentrate  in  a  small 
region  with  high  density.  When  the  beam  density  is  lower,  the  effect  is  minor,  an  influences  on  the  beam 
are  small.  For  higher  values  of  the  density,  the  effect  of  the  beam  is  much  stronger,  and  changes  in 
ambient  plasma  conditions  in  the  beam  region  has  a  severe  effect  on  beam  propagation.  One  can  probably 
consider  this  a  transient  stage,  where  a  channel  is  being  formed  in  the  plasma  allowing  later  portions  of 
the  beam  to  propagate  with  much  smaller  energy  loss. 
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Figure  3.7:  Beam  and  plasma  densities  on  the  .v-a.\is  as  a  function  of  time.  Left:  nb°lnp-5\  right 
«A7rtp=50.  Top:  beam;  bottom;  plasma. 


3.4  PIC  Summary 

Spacecraft  charging  effects  have  been  considered  looking  towards  experiences  with  keV  beams,  as  the 
principles  of  charging  are  not  expected  to  depend  on  beam  energy.  Charging  effects  are  expected  only  for 
ambient  densities  below  10’ m'^. 

Beam  dynamics  is  known  to  depend  critically  on  beam-to-ambient  plasma  densities,  and  thus  the  radial 
expansion  of  a  beam  is  of  primary  interest.  The  question  is,  will  a  beam,  injected  with  a  high  density, 
expand  to  low  densities  before  significant  beam-plasma  interactions  take  place?  Our  results  indicate,  that 
for  the  currents  considered  for  the  linac  beam,  there  will  be  no  problem  at  least  within  several  lO’s  of  km 
of  the  injection  point.  For  higher  current  beams,  where  ion-focusing  plays  a  role  (above  10  A),  beams 

will  attempt  to  form  an  ion  channel  in  the  plasma  and  will  expend  some  energy  in  this  process.  For  the 

38 


short  pulses  considered  here,  a  significant,  if  not  all  of  the  beam  pulse,  will  be  involved  in  this  process, 
while  for  msKiutation  beams,  a  laiger  portion  of  the  beam  may  propagate  relatively  stable  in  the  ion- 
channel  formed  by  the  beam  front. 

For  longer  distance  propagation,  the  beam  will  experience  scattering  in  both  energy  and  pitch  angle. 
Upward  injection  from  the  ionosphere  will  result  in  a  significant  portion  of  the  beam  being  scattered  out 
of  the  loss  cone,  to  form  curtains  of  energetic  electrons  drifting  around  the  Earth  like  radiation  belt 
electrons.  Estimates  of  energy  loss  in  the  process  is  of  the  order  of  10%  or  less,  a  number  found  for 
relativistic  electrons  injected  into  the  magnetosphere  from  above  thunderstorms  during  Sprite  generation 
[Lehtinen  et  al.,  2000].  For  beam  injection  downwards  into  the  upper  atmosphere  from  the  ionosphere,  a 
situation  of  interest  when  studying  properties  of  artificial  upward  discharges,  only  a  small  fraction  of  the 
beam  energy  will  be  lost  to  beam-plasma  interactions, 

4.  Time  Dependent  Relativistic  Electron  Beam  Distribution  Simulations  in 
the  Ionosphere-Magnetosphere;  Khazanov4D  and  KhazanovSD  Models 

This  task  was  accomplished  by  updating  time-dependent,  field-aligned  kinetic  models  of  electron 
transport  [Khazanov  et  al,  1994]  for  relativistic  energies  and  interactions.  Two  models  were  revised  and 
applied  to  the  present  situation.  The  first  is  a  4D  model  (one  space,  2  velocity  space,  and  one  time)  which 
we  refer  to  here  as  the  Khazanov4D  model.  A  second  model,  with  gyration-  and  bounce-averaging, 
referred  to  here  as  the  KhazanovSD  model  [Khazanov  et  al,  1996],  was  also  extended  to  relativistic 
energies  and  waves  were  included  in  the  model. 

The  following  summarizes  the  study  of  relativistic  electron  beam  injection  simulations  being  conducted 
with  the  Khazanov4D  and  Khazanov  5D  plasma  transport  models. 
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4.1.  Theory  and  Cross  Section  Discussion 

In  the  Kha2anov4D  model,  a  gyration-averaged  approach  is  assumed  because  the  electron  cyclotron 
period  is  by  far  the  smallest  time  scale  under  consideration.  In  this  limit,  the  field-aligned  kinetic 
equation  becomes 


vdt  ^  ds  IBdsdn^' 


(4.1) 


where  0  is  the  flux  of  electrons  in  units  of  cm’^s  'eV'sr  ',  E  is  the  kinetic  energy  of  the  electron  and 
fi=cos9.  Field-aligned  forces  have  also  been  omitted  from  this  equation,  because  it  is  expected  that  there 
will  be  no  parallel  forces  significant  enough  to  alter  the  relativistic  electron  trajectory.  The  left-hand  side 
of  (4.1)  represents  collisionless  transport  along  the  geomagnetic  field,  while  the  right-hand  side 
represents  all  collisional  processes  by  (S)  and  sources  by  Q.  Using  a  continuous  loss  approximation, 
the  atmospheric  electron-neutral  interactions  within  (5)  have  the  form 


d 


2  dpi 


(4.2) 


(4.3) 


Equation  (4.2)  represents  elastic  scattering  with  each  atmospheric  neutral  species  a,  where  is  the 
momentum  transfer  cross  section, 

=  J  4  X)(^  -  (4.4) 

where  /« is  the  Moliere  elastic  cross  section  discussed  in  previous  reports  and  %  is  the  scattering  angle. 
In  (4.3),  (7^  is  the  total  ionization  cross  section  for  species  a  into  ionization  state  r  calculated  by 
integrating  I^{E,E'),  the  Moller  inelastic  cross  section  discussed  in  previous  reports,  over  the 
secondary  electron  energies  E’,  E^  is  the  ionization  potential,  and  (E^)^  is  the  average  energy  of  the 

secondary  electron  spectrum. 
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(4.5) 


^Oar 


Also,  excitation  collisions  are  neglected  because  the  cross  sections  are  vastly  smaller  than  the  ionizaticr. 
cross  sections.  Therefore,  the  collisions  of  relativistic  electrons  with  atmospheric  neutral  particles  car. 
be  reduced  to  two  lookup  tables  for  and  +  (^s)ar)^^  ’  ^  functions  of  the 

r 

atmospheric  constituent  and  the  electron  energy.  The  first  quantity  is  the  momentum  transfer  cross 
section,  and  the  second  quantity  we  will  call  an  energy-loss  cross  section,  with  umts  of  eVcm  . 

We  will  begin  our  calculations  with  only  the  three  principal  atmospheric  species  included  in  the 
calculations:  O,  O2,  and  Nj.  The  momentum  transfer  cross  sections,  from  (4.4),  for  these  species  are 

shown  in  Figure  4.1. 

The  calculation  of  the  second  term  is  not  so  straight-forward,  because  information  must  be  knour. 
about  the  ionization  potentials  and  kinetic  energies  of  each  electron  in  the  atmospheric  constituen.:. 
First,  let  us  examine  the  dependence  of  the  cross  section  on  the  kinetic  energy  of  the  target  electro- 
(that  is,  the  electron  to  be  expelled  from  the  cloud  surrounding  the  nucleus). 

Figure  4r2  shows  results  for  a  binding  energy  of  13.6  eV  and  four  values  of  the  target  electron  kinetic 
energy:  0.01  eV,  1  eV,  13,6  eV,  and  150  eV.  As  seen  in  this  figure,  there  is  very  little  dependence  of  the 
cross  section  on  the  kinetic  energy  of  the  target  electron.  In  principal,  these  kinetic  energies  can  be 
calculated  theoretically  because  it  is  primarily  dependent  on  the  velocity  required  to  maintain  an  orbi: 
around  the  nucleus,  but  the  cross-electron  interactions  make  the  calculations  prohibitively  cumbersome, 
and  they  are  not  critical  to  the  result  in  any  case.  Therefore,  the  kinetic  energy  will  be  taken  as  equal  :c 
the  ionization  energy  for  that  target  electron. 

Uncertainty  is  still  found  in  the  choice  of  binding  energy,  however,  and  also  in  the  number  of  target 
electrons  associated  with  each  ionization  potential.  Test  cases  have  shown,  though,  that  this  is  muc.-. 
less  of  a  problem  than  expected  because  the  inner  shell  electrons  with  ionization  energies  more  than  an 
order  of  magnitude  larger  than  the  outer  shell  energies  have  cross  sections  far  below  the  outer  she;, 
results.  This  is  seen  in  Figure  4.3,  where  ionization  values  of  15  eV,  300  eV,  1000  eV,  and  1500  e\ 
have  been  used.  It  is  clear  that  only  the  ionization  of  the  outer  shell  \vill  significantly  contribute  to  the 
total  cross  section,  particularly  for  relativistic  energies,  even  when  the  inner  shell  cross  sections  are 
multiplied  by  the  number  of  target  electrons  in  these  shells. 
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Figure 4. 1:  Moliere  momentum  transfer  cross  sections  for  O  (solid  line),  O2  (dotted  line),  and  N2 
(dashed  line). 


C 

o 

a> 

{/) 

in 

(/> 

O 

w 

O 


Figure  4.2;  Inelastic  cross  sections  for  13.6  eV  ionization  energy  with  target  electron  kinetic  energies 
of  .01  eV  (solid  line),  1  eV  (dotted  line),  13.6  eV  (dashed  line),  and  150  eV  (dash-dot  line). 


It  is  also  clear  to  see  from  Figure  4.3  how  sensitive  the  cross  section  is  to  the  choice  of  the  ionization 
energy.  This  sensitivity  indicates  that  it  is  appropriate  to  take  into  account  the  various  excitation 
states  in  which  the  newly  created  ion  can  be  left.  Values  for  these  ionization  potentials,  and  the 
corresponding  branching  ratios  for  each,  are  given  in  Table  4.1.  These  values  are  from  the  non- 
relativistic  ionization  potentials  given  by  Jackman  et  al.  [1977].  The  branching  ratios  were  determined 
from  the  high-energy  range  of  the  validity  of  these  cross  sections,  and  are  multiplied  by  the  number  of 
electrons  in  the  outer  shell  of  each  constituent,  taken  as  2, 4,  and  4  for  the  three  species  of  interest  (0, 
O2,  and  N2,  respectively). 
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Figure  4.3:  Inelastic  cross  sections  with  ionization  energies  of  15  eV  (solid  line),  300  eV  (dotted  line), 
1000  eV  (dashed  line),  and  1500  eV  (dash-dot  line). 


Table  4.1 

;  Ionization  states  and  branching  ratios 

0  states 

E" 

Ratio 

O2 

states 

Ratio 

N2 

states 

E"" 

Ratio 

■‘s 

13.60 

.370 

X 

12.10 

.079 

X 

15.58 

.456 

16.90 

.410 

a 

16.10 

.191 

A 

16.73 

.200 

‘P 

18.50 

.210 

A 

16.90 

.191 

B 

18.75 

.104 

b 

18.20 

.170 

D 

22.00 

.070 

B 

20.00 

.110 

C 

23.60 

.070 

c 

23.00 

.160 

40 

40.00 

.100 

37 

37.00 

.099 

Cross  section  results,  along  with  the  average  seconda^  electron  energy  from  (4.5)  and  the  total 
"energy-loss  cross  section,”  are  presented  in  Figures  4v4-4.6.  In  the  cross  section  plots  (top  panels), 
the  solid  line  is  the  summation  of  the  partial  cross  sections  for  each  ionization  state.  Note  that  these 
cross  sections  are  integrated  over  the  range  of  secondary  electrons,  and  that  all  three  species  have  a 
cross  section  of  about  2x10'^°  cm^  for  relativistic  primary  energies.  The  middle  panels  show  the 
average  secondary  electron  energies  for  each  ionization  state,  and  the  solid  line  is  a  weighted  average, 

^^ar{^e)ar  (4.6) 
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It  is  interesting  that  the  average  secondary  energy  from  a  relativistic  primary  energy  is  about  300  eV, 
while  the  average  secondary  in  the  non-relativistic  limit  is  about  20  eV.  In  the  lower  panel,  the  solid 
line  is  the  result  from  summing  the  various  ionization  states  while  the  dotted  line  is  the  result  with  all 
of  the  newly  created  ions  starting  in  the  ground  state.  These  results  show  that  assuming  all  of  the  outer 
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shell  electron  ionizations  jump  to  the  ground  state  of  the  ionized  particle  is  valid,  and  that  using  this 
table  of  branching  ratios  into  the  various  ionization  states  is  merely  fine-tuning  to  the  overall  quantity 
needed  in  the  model. 


Figure  4.4:  Cross  sections,  average  secondary  energies,  and  total  energy-loss  cross  section  for  atomic 
oxygen. 
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Figure  4.6:  Cross  sections,  average  secondary  energies,  and  total  energy-loss  cross  section  for 
molecular  oxygen. 


4.2  Khazanov4D  Model  Simulations 

The  Kha2anov4D  model  was  modified  to  include  relativistic  effects.  Several  test  cases  were  conducted 
to  ensure  proper  simulation,  and  an  example  case  is  shown  below.  Results  using  the  model  are  also 
reported  in  Khazanov  et  al.  [1999] 

For  the  results  below,  the  model  was  run  with  one  energy  step  0.25  MeV  wide  centered  at  5  MeV, 
injected  at  700  km  at  pitch  angles  greater  than  40”  (upward  direction  only).  The  injection  was  assumed 
to  be  Gaussian  with  a  width  of  0.5  MeV.  It  was  run  with  a  time  step  of  lO’^  s  for  an  L=2  field  line 
(ionospheric  latitude  near  45”,  with  an  equatorial  plane  crossing  at  2  Rg  geocentric).  The  injection 
lasted  one  time  step,  and  then  the  beam  was  allowed  to  decay  with  time.  This  configuration  was 
chosen  because  it  illustrates  the  longest  lifetime  of  a  relativistic  beam  injection  fi-om  the  ionosphere. 
That  is,  by  injecting  at  700  km,  particles  with  large  pitch  angles  will  mirror,  or  bounce  (due  to  the 
inhomogeneous  pomagnetic  field),  in  the  upper  ionosphere,  before  they  experience  excessive  collisions 
with  atmospheric  neutral  particles  (which  we  feel  will  be  the  largest  loss  mechanism  fdr  these 
particles).  Also,  because  the  simulations  do  not  include  the  effects  of  plasma  waves  (either  externally- 
generated  waves,  like  hiss  and  whistlers,  or  self-generated  waves  due  to  the  unstable  nature  of  the  beam 
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(4.8) 
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The  solid  line  is  the  same  run  as  shown  in  the  previous  plots,  while  the  dotted  line  shows  the  results 
from  a  run  with  no  scattering  or  loss  mechanisms  included  (except  loss  to  altitudes  below  200  km).  It 
is  seen  that  after  the  initial  loss  of  particles  due  to  the  injection  pitch  angle  range  spanning  part  of  the 
loss  cone,  the  model  does  not  have  significant  loss  of  particles.  Also  seen  is  the  exponential  decrease  in 
the  number  of  particles  as  seen  in  Figures  4.7  and  4.8.  The  half-life  of  the  beam  particle  count  is  0.68  s. 
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Figure  4.7:  Downward  fluxes  at  the  injection  location  (700  km  altitude  in  the  northern  ionosphere)  as 
the  beam  passes  by  this  spatial  point.  The  results  are  for  successive  bounces  of  the 
beam,  with  the  heavy  line  being  the  initial  distribution  and  a  solid  line  every  5th  bounce. 
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Figure  4.10:  Equatorial  pitch  angle  distributions  during  successive  bounces  of  the  injected  beam.  The 
upper  panel  shows  results  for  the  first  2  s  with  a  curve  every  bounce  period,  and  the 
lower  panel  shows  results  for  the  next  8  seconds,  with  a  curve  every  5th  bounce.  The 
trapped  zone  flux  is  increasing  in  the  upper  panel,  but  is  then  decreasing  in  lower  panel. 


Figure  4.11:  Total  beam  particle  evolution  for  a  run  with  scattering/losses  (solid  line)  and  without 
scattering/losses  (dotted  line). 

4.3  KhazanovSD  Injection  Configuration 

The  configuration  of  the  injected  relativistic  electron  beam  is  the  same  as  that  used  for  the  Khazanov4D 
model,  as  discussed  in  Khazanov  et  al  [1999].  It  assumes  an  upward  injection  from  700  km  altitude  in 
the  ionosphere,  with  particle  loss  to  the  atmosphere  occurring  at  200  km.  The  beam  is  injected  with  a 
Gaussian  energy  spectra  centered  on  5  MeV  with  a  0.5  MeV  width,  and  an  intensity  that  is  constant  with 
pitch  angle.  The  energy  grid  is  divided  into  9  cells  0.25  MeV  wide,  the  lowest  one  centered  on  4  MeV 
and  the  highest  one  centered  on  6  MeV.  The  equatorial  pitch  angle  grid  is  chosen  so  that  5  grid  cells  fit 
exactly  within  trapped  part  the  beam  injection  region  (that  is,  that  part  of  the  beam  that  mirrors  above  200 
km),  and  then  the  cell  size  smoothly  enlarges  toward  0"  and  90'.  For  L=2,  3,  and  4,  this  trapped 
injection  range  is  2.33',  1.13',  and  0.67'  wide  centered  about  18.37°,  9.39',  and  5.95',  respectively. 
For  the  initial  MLT  distribution,  a  Gaussian  centered  about  local  noon  with  a  0.25  hour  spread  was 
assumed.  This  is  not  only  for  numerical  stability,  avoiding  sharp  gradients  in  the  variable  of  greatest 
advective  drift,  but  also  to  simulate  the  initial  spread  of  the  injected  beam  over  1  bounce. 
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43.1  Model  Parameters 

The  collisional  operators  in  the  kinetic  equation  require  information  about  the  neutral  atmosphere  and 
thermal  plasma.  For  these  background  quantities,  the  standard  models  of  MSIS-90  and  IRI-90  are  used 
at  spring  equinox  of  1991.  This  represents  a  solar  maximum  condition  when  collisional  losses  are 
expected  to  be  greatest.  The  equatorial  plane  thermal  plasma  densities  are  taken  from  a  quiet  time  result 
from  the  Rasmussen  plasmaspheric  model,  which  is  the  standard  dynamic  plasmasphere  model  used  in 
the  KhazanovSD  model. 

Wave-particle  interaction  diffusion  coefficients  are  calculated  using  the  method  of  Albert  [1999].  Four 
wave  types  are  included  in  the  wave-particle  diffusion  coefficient  calculation:  plasmaspheric  hiss; 
lightning-generated  whistlers;  and  two  bands  of  VIP  transmitter  signals.  All  are  right-hand  polarized 
whistler-mode  waves  with  Gaussian  distributions  in  both  frequency  and  wave  propagation  angle.  The 
characteristics  used  for  these  waves  are  those  given  by  Abel  and  Thome  [1998]. 

These  background  particle  and  wave  characteristics  are  used  to  calculate  pitch  angle  scattering  and  energy 
loss  coefficients.  A  representative  sample  of  the  diffusion  coefficients  is  shown  in  Figure  4.12  at  several 
spatial  locations.  The  scattering  rates  for  each  type  of  wave  as  well  as  that  for  particle  collisions  are 
presented,  as  well  as  the  initial  condition  pitch  angle  range  for  reference.  It  is  clear  that  collisional 
scattering  is  only  effective  for  pitch  angles  that  reach  the  atmosphere,  with  this  coefficient  rising  by  two 
orders  of  magnitude  over  its  plasmaspheric  value.  The  wave  coefficients  are  dominated  everywhere  by 
the  interaction  with  hiss,  and  drop  significantly  at  lower  pitch  angles.  However,  the  scattering  rate  at  the 
inner  edge  of  the  injection  range  is  roughly  equal  for  the  two  processes  at  L=2  and  3,  and  they  are  within 
a  factor  of  2  at  L=4.  The  roughness  in  the  VLF  coefficients  is  due  to  the  narrowness  of  the  frequency 
bandwidth  for  these  waves. 
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Figure  4.12:  Magnitudes  of  the  bounce-averaged  pitch  angle  diffusion  coefiBcients  versus  equatorial 
pitch  angle  at  1=2,  3,  and  4  and  at  MLT=midnight  and  noon  for  5  MeV  electrons 
interacting  with  hiss  (solid  line),  lightning  whistlers  (dotted  line),  the  VLFl  band  (dashed 
line),  the  VLF2  band  (dash-dot  line),  all  waves  (dash-dot-dot-dot  line),  and  particles  (long 
dashed  line).  The  small  black  boxes  indicate  the  piteh  angle  injection  range. 

A  sample  of  the  collisional  energy  loss  rates  is  given  in  Figure  4vt3  at  several  spatial  locations  and 
several  pitch  angles  near  the  injection  regions.  Again  note  the  dramatic  increase  (of  2+  orders  of 
magnitude)  in  loss  rate  for  pitch  angles  mirroring  in  the  atmosphere.  The  dashed  line  is  at  the  center  of 
the  initial  injection  range,  yet  the  loss  rate  is  minor  compared  to  that  at  the  loss  cone  edge.  Figures  4.-13 
and  4.14  will  be  discussed  in  more  detail  below. 
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Figure  4.13:  Magnitudes  of  the  bounce-averaged  energy  loss  coefficients  for  several  equatorial  pitch 
angles  at  Z,=2,  3,  and  4  and  at  MLT=noon  and  midnight. 


4.3.2  Global  Representation 

In  order  to  examine  the  global  evolution  Of  an  injected  relativistic  electron  beam  in  the  inner 
magnetosphere,  simulations  were  conducted  using  the  above-mentioned  initial  and  background 
conditions  for  L=2,  3,  and  4  (41.8°,  52.5°,  and  58.2°  magnetic  latitude  at  700  km).  The  gyration- 
averaged,  bounce-averaged,  time-dependent,  relativistic  kinetic  equation  was  solved  for  the  distribution 
function  of  the  beam  electrons  as  a  function  of  time,  azimuth  (MLT),  energy,  and  equatorial  pitch  angle. 
Because  the  calculation  is  linear,  the  results  presented  below  are  given  in  normalized  units  and  can  be 
scaled  to  match  any  beam  injection  intensity. 


Figure  4.14  shows  the  distribution  function  versus  energy  and  MLT  for  the  pitch  angle  cell  at  the 

midpoint  of  the  initial  condition  region  for  L=2.  These  results  are  for  simulations  with  all  processes 

included  (drift,  collisions,  and  wave  interactions).  The  contours  are  spaced  every  half  order  of 

magnitude,  and  the  normalization  is  scaled  to  the  maximum  value  at  t=0.  As  the  simulation  progresses, 

the  beam  splits  into  a  series  of  packets  (at  all  three  L  values).  This  is  a  numerical  effect  from  choosing 
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only  9  energy  grid  cells.  The  dotted  line  through  the  packets  demarks  the  drift  pattern  at  each  time.  With 
infinitely  many  energy  grid  cells,  the  result  at  t>0  would  look  similar  to  the  result  at  ts=0  but  tilted  along 
this  dotted  line  and  collisionally  spread.  After  two  drift  periods,  the  6  MeV  electrons  are  about  to  lap  the 
4  MeV  electrons,  and  a  banded  structure  appears  in  the  energy  spectra.  While  the  simulation  results  still 
have  a  beaded-string  appearance  along  the  dotted  line,  the  reality  (replicated  by  interpolating  between  the 
beads  along  the  dotted  line)  is  a  multiply-spiked  energy  spectra  at  all  local  times  (banded  torus).  After  a 
few  hours,  the  overlap  is  severe.  The  packets  are  also  spreading  due  to  collisions  and  scattering.  At  6 
hours,  the  packets  overlap  to  such  a  degree  that  interpolating  them  along  the  dotted  lines  (as  the  beam  is 
in  reality)  forms  a  non-banded  torus  around  the  Earth.  As  the  simulation  progresses,  the  packets 
eventually  disappear  (even  in  spite  of  the  coarse  energy  grid)  and  the  beam  has  become  (numerically  and 
physically)  a  smooth  torus  of  relativistic  electrons  at  the  L  shell  of  injection. 


Figure  4. 14:  The  beam  distribution  versus  energy  and  MLT  for  1 8.4'  and  L=1  at  various  times. 

Figure  4.15  shows  the  pitch  angle  versus  MLT  distributions  at  5  NfeV  for  the  same  simulations  as  in 
Figure  4.14.  These  results  are  also  normalized  to  the  peak  value  of  the  initial  condition  (t=0).  The  initial 
condition  shows  a  narrow  Gaussian  in  MLT  and  uniform  distribution  across  the  pitch  angle  range 
between  the  pitch  angles  that  mirror  between  200  km  in  one  conjugate  hemisphere  and  700  km  in  the 
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opposite  hemisphere.  Because  the  drift  differential  across  the  pitch  angle  injection  range  is  less  than  for 
the  energy  range,  the  beam  does  not  split  into  packets,  and  actually  appears  as  a  tilted  version  of  the 
initial  distribution.  An  interesting  feature  develops  in  addition  to  this  tilting  effect:  a  leading  population 
forms  right  at  the  edge  of  the  200  km  loss  cone.  This  feature  lengthens,  broadens,  and  intensifies  with 
time.  This  feature  is  due  to  collisional  energy  losses.  This  term  is  strongest  right  near  the  edge  of  the 
loss  cone  where  the  electrons  mirror  deep  in  the  atmosphere  and  encounter  a  laige  number  of  scattering 
targets.  As  the  higher  energy  electrons  continue  to  outpace  the  5  MeV  electrons  shown,  they  continue  to 
deposit  decayed  electrons  at  local  times  farther  and  farther  ahead  of  the  beam  location.  The  beam, 
meanwhile,  is  diffusing  in  pitch  angle  (and  thus  in  MLT  via  the  drift  differential),  eventually  becoming  a 
broad  band  of  electrons  at  all  MLT.  After  one  day,  this  band  is  more  than  100%  wider  than  the  initial 
injection  width.  Because  of  the  smallness  of  the  pitch  angle  drift  differential,  however,  significant 
structure  still  exists  in  this  band,  with  an  oblique  ridge  as  a  distinctive  remnant  of  the  original  injection 
configuration. 


Figure  4. 15:  The  beam  distribution  versus  equatorial  pitch  angle  and  MLT  for  5  MeV  and  Z,=2  at 
various  times. 
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4.33  Total  Beam  Particle  Content 

With  so  much  differential  spreading  over  MLT,  it  is  difficult  to  extract  the  magnitude  of  the  beam  loss  as 
a  function  of  time  from  plots  of  the  distribution  function.  A  useful  quantity  to  examine,  therefore,  is  the 
total  number  of  particles  in  the  beam.  This  quantity  is  found  by  integrating  the  distribution  function  over 
all  phase  space  at  a  given  time.  Figure  4.16  presents  this  quantity  for  various  simulations.  The  upper 
panel  shows  results  at  L=2  for  azimuthal  transport  (drift)  only,  drift  plus  collisions,  drift  plus  wave 
interactions,  and  all  processes  included.  First,  the  simulation  for  drift  only  solidly  maintains  its  original 
number  of  particles.  This  is  encouraging,  because  atmospheric  loss  was  not  omitted  from  this  calculation 
(it  is  included  in  all  of  these  simulations),  and  the  scattering  and  loss  operators  were  also  not  omitted, 
simply  calculated  with  coefficients  of  zero.  This  means  that  there  was  no  numerical  diffusion  from  these 
calculations.  This  is  also  seen  in  plots  similar  to  Figures  4.14  and  4.15  (not  plotted,  however),  which 
show  the  simulation  maintaining  the  original  functional  form  of  the  beam.  Secondly,  this  plot  shows  that 
waves  are  more  important  at  this  L  shell  than  collisions  with  particles  for  removing  electrons  from  the 
beam,  but  collisions  are  not  an  insignificant  process.  Note  that  it  is  not  a  linear  addition  of  the  two  loss 
amounts  to  yield  the  result  when  all  processes  are  included,  but  there  is  actually  less  loss  in  this  case. 
This  is  because  there  is  only  a  limited  number  of  particles  in  the  initial  beam,  and  once  they  are  scattered 
into  the  loss  cone,  no  matter  what  process  does  it,  they  are  gone  from  the  simulation. 

The  lower  panel  of  Figure  4.16  shows  the  total  particle  count  for  simulations  with  all  processes  included 
at  various  L.  While  the  beam  loses  35%  of  its  particles  after  24  hours  at  L=2,  it  only  loses  30%  at  L=3, 
but  loses  60%  at  L=4.  The  decay  rates  at  t=24  h  for  the  three  simulations  are  0.99%/h,  0.94%/h,  and 
1.75%/h,  which  are  all  slowly  decreasing  with  time.  Compare  these  late-stage  loss  rates  to  the  loss  rates 
1  hour  after  injection:  2.4%/h,  2.6%/h,  and  13.8%/h  for  the  3  L  values,  respectively.  These  results  can 
be  explained  by  reexamining  Figures  4.12  and  4.T3.  The  coefficients  in  those  plots  can  be  roughly 
converted  into  time  scales  of  particle  energy  loss  or  scattering  through  an  £  or  interval.  Right  at  the 
loss  cone  edge,  it  should  take  about  0.8,  1,  and  1.4  hours  for  an  electron  to  lose  1  MeV  of  energy  for 
L=2,  3,  and  4,  respectively  (note  that  it  takes  orders  of  magnitude  longer  at  pitch  angles  away  fi’om  the 
loss  cone  edge).  The  scattering  time  scale  at  the  loss  cone  edge  is  about  3.5, 3.3,  and  0.6  hours  for  L=2, 
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3,  and  4,  respectively.  The  dramatic  drop  for  L=A  is  because  the  width  of  the  pitch  angle  injection 
shrinks  with  L  faster  than  the  diffusion  coefficient.  These  time  scales  are  crude  estimations  for  a  single 
energy  and  pitch  angle,  and  not  global  decay  rates  for  all  phase  space. 


The  lower  panel  of  Figure  4.16  can  now  be  viewed  in  terms  of  the  relative  magnitudes  of  these  loss 
rates.  The  initial  loss  rate  is  highly  dependent  on  the  scattering  rate  at  the  edge  of  the  loss  cone,  because 
this  is  a  sharp  gradient  next  to  a  precipitation  zone.  The  order  of  the  initial  loss  rates  of  the  L  shells  match 
our  scattering  time  scale  estimations.  The  late-stage  loss  rate  is  caused  by  a  combination  of  scattering  and 
decay  losses.  While  L=3  has  a  similar  diffusion  loss  rate  to  that  of  L=2,  the  energy  loss  rate  is  smaller 
and  thus  the  beam  retains  more  of  its  particles.  1=4  however  has  a  much  faster  diffusive  loss  rate  that 
overcomes  the  slower  energy  loss  rate  to  result  in  a  larger  late-stage  loss  rate  (and  reduced  beam 
intensity)  at  this  L  value. 


Figure  4.16:  Evolution  of  the  normalized  total  number  of  particles  of  the  beam  at  L=2  for  several 
simulations  with  various  processes  included  (upper  panel)  and  at  several  L  values  with  all 
processes  included  (lower  panel). 
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4.4  Khazanov  Modeling  Summaiy 

In  order  to  examine  the  global  evolution  of  an  artificially  injected  relativistic  electron  beam,  the  bounce- 
averaged  relativistic  kinetic  equation  was  solved  for  the  electron  distribution  function  for  various 
scenarios.  It  was  found  that  the  beam  quickly  spreads  in  MLT  due  to  differential  drift  rates.  This, 
combined  with  energy  loss  and  pitch  angle  diffusion,  eventually  transform  the  beam  into  a  fairly  uniform 
torus  around  the  Earth.  It  was  determined  that  wave  interactions  are  slightly  more  important  than 
collisional  losses  in  reducing  the  1=2  beam  content.  It  was  also  found  that  the  beam  total  particle  loss 
rate  is  complicated  function  of  L  and  initial  condition  parameters,  and  for  the  chosen  conditions,  a  beam 
at  1=4  drops  much  faster  than  a  injected  beam  at  L=2  or  3.  These  loss  rates  are  about  1-2%/h  one  day 
after  injection,  with  the  beam  total  particle  counts  at  65%,  70%,  and  40%  of  the  initial  count  at  this  time. 
Other  interesting  features  in  the  distribution  function  were  examined,  such  as  the  leading  edge  population 
in  the  pitch  angle-MLT  plots,  banded  structures  in  the  energy  spectra,  and  reasons  for  the  values  of  the 
early  and  late  stage  decay  rates. 

A  final  note  of  discussion  is  on  the  comparison  of  total  loss  rates  determined  in  this  study  with  those  of 
our  previous  study  [Khazanov  et  al.,  1999].  There,  a  time  scale  of  265  s  was  stated  for  an  e-folding  of 
the  beam  intensity.  However,  because  of  the  very  small  time  step  required  to  resolve  the  field-aligned 
structure  of  the  beam,  those  simulations  were  only  conducted  for  the  first  10  s  after  injection.  As 
indicated  above,  the  loss  rate  is  initially  quite  high  but  then  dramatically  decreases  during  the  first  few 
hours  after  injection.  In  fact,  the  decay  rates  for  all  simulations  are  still  declining  after  24  hours.  Thus, 
the  Khazanov  et  al  [1999]  value  is  a  good  estimate  of  the  initial  decay  rate  of  the  beam,  while  this  study 
has  examined  a  different  time  regime. 

Later  studies,  based  on  this  relativistic  model  development,  have  also  shown  that  the  beam  particles 
strongly  depend  on  the  initial  pitch  angle  distribution  and  that  lifetimes  similar  to  those  found  for  the 
radiation  belts  are  obtained  for  nearly  equatorial  minroring  injections  [Khazanov  et  al.,  20(X)].  It  was 
concluded  that  the  Khazanov  et  al.  [1999]  results,  which  only  considered  upper  ionospheric  injections, 
are  consistent  with  other  relativistic  electron  studies  and  should  be  regarded  as  an  examination  of  the  loss 
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cone  edge  population,  whether  naturally  occurring  or  artificially  injected  into  near-Earth  space.  It  is  also 
noted  that  the  model  is  well  adapted  to  the  analyses  of  the  radiation  belt  fine  structure  evolution  and  the 
dynamics  of  various  regions  of  velocity  space. 


5.  Summary,  Recommendations,  and  Publications 

5.1  Summary  of  Accomplishments  and  Results 

The  following  represent  key  accomplishments  and  results  for  this  project: 


•  The  RelBAI  model  was  developed  to  allow  detailed  accounting  of  beam  energy  propagation, 
generation  of  secondary  particles,  and  radiation  effects  such  as  bremsstrahlung; 

•  The  RelBAI  model  was  used  to  complete  analysis  of  representative  cases  associated  with 
propagating  a  relativistic  electron  beam  from  low-Earth  orbit  down  in  to  the  atmosphere; 

•  Dr.  Linda  Habash  Krause  completed  her  Ph.D.  based  on  the  development  and  application  of  the 
RelBAI  model; 

•  The  RelBAI  model  was  extended  beyond  the  project  work  described  here  to  include  a  relativistic 
electron  auroral  model  to  address  the  time-dependent  chemistry  necessary  to  compute  optical 
emissions  as  a  function  of  altitude  and  several  potential  emissions  were  identified; 

•  3D  PIC  simulation  models  were  developed  for  the  relativistic  beam  case; 

•  The  PIC  modeling  indicates  that  charging  effects  will  only  be  an  issue  at  ambient  plasma  densities 
less  than  lO’ m*^  for  the  proposed  beam  current  and  duty  cycle; 

•  Ion  focusing  of  the  beam  will  play  a  role  only  for  beam  currents  above  10  A,  which  is  well  above 
the  proposed  beam  current; 

•  PIC  results  indicate,  that  for  the  currents  considered  for  the  linac  beam,  there  will  be  no  significant 
beam-plasma  interaction  creating  undesirable  beam  dynamics  at  least  within  several  IQ’s  of 
kilometers  of  the  injection  point; 

•  For  longer  propagation,  the  PIC  results  suggest  the  beam  wiU  e>q)erience  scattering  in  both  energy 
and  pitch  angle  though  it  will  continue  to  propagate; 

•  Khazanov4D  and  KhazanovSD  models  were  extended  to  relativistic  beam  energies; 

•  The  Khazanov  models  showed  that  the  relativistic  electron  beam  quickly  spreads  in  MLT  due  to 
differential  drift  rates  and  combined  with  energy  loss  and  pitch  angle  diffusion,  eventually 
transform  the  beam  into  a  fairly  uniform  torus  around  the  Earth; 

•  The  Khazanov  models  showed  that  the  beam  total  particle  loss  rate  is  complicated  function  of  L 
and  initial  condition  parameters; 
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•  The  Khazanov  models  showed  that  the  beam  particles  strongly  depend  on  the  initial  pitch  anglp 
distribution  and  that  lifetimes  similar  to  those  found  for  the  radiation  belts  are  obtained  for  nearly 
equatorial  mirroring  injections. 

5.2  Recommendations  for  Future  Work 

It  is  recommended  to  go  forw^d  wiSi  the  relbe^  flight  progranu  As  a  first  step,  it  is  proposed  to 
develop  an  accelerator  fli^t  unit  for  a  sounding  rocket  mission.  During  flight,  the  accelerator  should 
inject  beam  pulses  downwards  into  the  upper  atmosphere.  Simultaneous|^rdinated  oblations’  fioin 
ground  should  be  done  of  optical  emissions  and  radar  echoes  firom  density  striations.  This  scenario 
provides  for  the  best  diagnostics  configuration  of  accelerator  performance  and  beam  propagation 
characteristics. 

It  is  further  recommended  to- carry  out  the  experiment  in  connection  with  severe  thunderstorms,  in 
order  to  explore  electrical  processes  in  connection  \vith  thunderstorm  electric  charge-  and  discharging. 
The  phenomenon  of  sprites  is  now  understood  to  the  point  where  good  predictions  can  be  of 
time  and  location  of  sprites.  Typically,  sprites  occur  from  well  defined  thunderstorm  cells  over  a 
period  of  30  minutes  or  more,  during  which  time  the  occurrence  rate  is  about  one  sprite  every  2 
minutes.  The  details  of  sprite  generation  is  not  know,  out  it  is  thought  that  a  discharge  is  initiated  at  70 
km  altitude  in  connection  with  a  positive  cloud-to-ground  discharge  in  the  storm  below.  Possibly,  the 
discharge  process  can  be  explored  with  relativistic  beams,  where  seed  ionization  and  electric  potential 
modification  is  generated  down  to  40-50  km  altitude.  It  is  possible,^  a  relbeam  will  initiate  a  discharge. 
Such  experiment  would  be  unique  and  quite  spectacular,  if  successful.  It  would  represent  an  addfyl 
science  value  to  the  basic  experiment  objectives  that  are  to  assess  accelerator  and  beam  performance  in 
space. 

Launch  cites  must  be  assessed  for  their  compatibility  with  the  sprite  experiments.  It  is  required,  that 
the  beam  track  at  50  km  altitude  is  across  a  regions  with  severe  thunderstorm  activity.  White  Sands, 
New  Mexico,  and  Wallops  Island  may  be  candidates. 

The  theoretical  results  presented  here  have  shown  from  various  points  of  view  that  relativistic  beam 
propagation  in  the  ionosphere  and  magnetosphere  appears  feasible.  By  this  we  mean  that  the  beam  is 
expected  to  easily  escape  a  spacecraft,  that  the  beam  appears  generally  stable  as  it  propagates  initially, 
and  it  is  possible  for  the  beam  to  have  very  large  lifetimes.  If  the  beam  pitch  angle  is  low  enough  to  be 
within  the  loss-cone,  the  beam  can  be  deposited  in  the  atmosphere. 

The  models  developed  under  this  project  are  ready  and  able  to  focus  in  on  specific  potential 
applications.  They  can  be  used  to  address  practical  aspects  or  issues  of  conducting  actual  flight 
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experiments.  Concerns  with  detection  methods  on  the  one  hand  and  impacts  to  other  space  assets  can 
both  be  addressed.  This  would  be  the  next  logical  steps  in  this  work. 


5.3  Publications  and  Presentations 

The  following  publications  and  presentations  were  generated  as  a  result  of  the  work  sponsored  by  this 
grant. 

PUBLICATIONS 

Khazanov,  G.  V.,  M.  W.  Liemohn,  E.  N.  Krivorutsky,  J.  U.  Kozyra,  and  B.  E.  Gilchrist, 
Interhemispheric  transport  of  relativistic  electron  beams,  Geophys.  Res.  Lett.,  26, 581, 1999. 
Neubert,  T.  and  B.  E.  Gilchrist,  3D  electromagnetic  PIC  simulations  of  relativistic  electron  pulse 
injections  from  spacecraft,  Adv.  Space  Res.,  submitted,  2000. 

Khazanov,  G.  V.,  M.  W.  Liemohn,  E.  N.  Krivorutsky,  J.  U.  Kozyra,  J.  M.  Albert,  and  B.  E. 
Gilchrist,  On  the  influence  of  the  initial  pitch  angle  distribution  on  relativistic  electron  beam 
dynamics,  J.  Geophys.  Res.,  105,  16093,  2000. 


PRESENTATIONS 

Habash  Krause,  L.,  Analysis  of  Active  Space  Experiments  Using  Artificial  Relativistic  Electron  Beams 

Habash  Krause,  L.,  T.  Neubert,  and  B.  E.  Gilchrist,  Passive  remote  sensing  of  artificial  relativistic 
electron  beams  in  the  middle  atmosphere,  Proc.  AIAA  1999  Space  Technology  Conference  and 
Exposition,  Norfolk,  VA,  AIAA-99-4352,  1 999. 

Habash  Krause,  L.,  B.  E.  Gilchrist,  and  T.  Neubert,  Analysis  of  active  space  experiments  with 
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